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We consider the propagation of hydromagnetic waves generated by a compact turbulent source in low 
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ionized plasmas, applying the Lighthill theory. We assume the plasma to be isothermal, and adopt a 
uniform, stationary medium thread by ordered magnetic fields as an initial condition. Then, the distinct 
properties of the hydromagnetic waves originating from a source oscillating with a fixed frequency are 
studied in the linear regime. As is well known, in low ionized plasmas, the generated waves dissipate due to 



ion-neutral damping. In this paper, the dependence of the dissipation rate on the frequency of the oscillating 
source is investigated. The larger the frequency becomes, the more substantial is the wave dissipation. 
Implications of our results on the energy source in molecular clouds are also discussed. Interestingly, since 
the outflow lobes associated with young stellar objects act as compact turbulent sources, hydromagnetic 
waves are generated by them. From our order-estimations, about 70% of the energy of the outflow itself 
propagates as waves or turbulences, while the remaining 30% dissipates and heats the neutrals via ion- 
neutral damping. Then, we confirm that the outflows are significant energy sources in molecular clouds in 
the context of the Lighthill theory. 
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Young stellar objects (YSOs) are believed to be formed within dense clumps in molecular clouds (Myers, Benson 
1983). Going with the gas condensation associated with star formation, a large amount of gravitational energy must 
be released, and a part of the released energy drives the jets and/or outflows (e.g., Shu et al. 1993; Kudoh, Shibata 
1997; Tomisaka 1998). In this paper, we mainly discuss the role of the outflows. The total mechanical luminosities 
of all the outflows in a molecular cloud are observed, for example, to be ~ 33 L m in Orion A (Fukui et al. 1993) 
and ~ 43.5 in Taurus (Bontcmps et al. 1997). We stress the importance of the activity of the outflows on the 
ambient molecular gas. 

Lighthill (1952) discussed the aerodynamical generation of waves from a compact turbulent source; his theory 
was developed concerning the generation of hydromagnctic waves in a magnetized medium (Kulsrud 1955; Kato 
1967). Waves originate from the pressure fluctuations caused by fluctuations in the momentum of velocity and 
magnetic fields of unsteady turbulent flows. Their theories have been applied to studies of the heating mechanism of 
stellar chromosphere and coronae (see Lee 1993 and references therein). In this paper, we apply them to a study of 
molecular clouds, where the ionization fraction is low. In star-forming regions, the outflow lobe filled with turbulent 
molecular gas (e.g., Stahler 1994) is one of the good candidates for a compact turbulent source (e.g., Nomura et 
al. 1998). We then examine the potentiality of outflows as an energy source in molecular clouds in the context of 
the Lighthill theory. Part of the energy of the outflows dissipates in the molecular gas via the so-called ion-neutral 
damping (Kulsrud, Pearce 1969; Zweibel, Josafatsson 1983; Nakano 1998, see also Kamaya, Nishi 1998), as far as 
the outflows can generate the hydromagnetic waves. 

The problem concerning the heating source in molecular clouds is of interest. In molecular clouds (with a temper- 
ature of ~ 10 K and a number density of ~ 10 3 cm -3 ) the cooling processes have been well studied (e.g., Goldsmith, 
Langer 1978; Neufeld et al. 1995). The heating processes, however, are problematic because it is uncertain what is 
the dominant process. Although it is often assumed that the main heating source is cosmic rays, the heating rate 
due to H 2 formation on grains, gas-grain collisions, compressions or collapses of dense cores, and turbulence may 
be comparable to that caused by cosmic rays (Lizano, Shu 1987; Hollenbach 1988; Stone et al. 1998). In this paper, 
we insist that the ion-neutral damping of waves generated from outflows is also a dominant heating mechanism (cf. 
Zweibel, McKee 1995). 
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We also discuss the role of outflows as an energy source to support the turbulence in molecular clouds. In many 

molecular clouds and cloud cores, the scaling law of the line width-size relation, Av oc L 1 / 2 , has been observed with 

various molecular lines (Larson 1981; Falgarone et al. 1992; Caselli, Myers 1995). The observed widths are wider 

than the thermal width and the power-law index of the scaling law is similar to that of the Kolmogorov turbulence. 

These facts suggest that the interstellar medium may be turbulent (Larson 1981; Scalo 1987; Franco, Carraminana 

1998). If so, it should quickly dissipate because of the shock, etc. (e.g., Stone et al. 1998); thus, some kinds of energy 

sources are needed. Many candidates have been suggested, e.g., galactic differential rotations, isolated supernovae, 

stellar winds, Alfven waves from collapsing cloud cores, expanding H II regions, and bipolar outflows (Norman, Silk 

1980; Mouschovias, Paleologou 1980; Fleck 1981; Scalo 1987; Miesch, Bally 1994). In recent numerical simulations it 

is suggested that the sources with stellar sizes may be more important than larger-scale sources (Gammie, Ostrikcr 

1996; Vazquez-Semadeni et al. 1997). In the following discussions, we estimate how much energy is given to the 

turbulence through waves generated by outflows associated with YSOs. 

In subsequent sections, we investigate the evolution of hydromagnetic waves generated from a compact turbulent 

source in a low-ionized, magnetized medium. We note that the difference between ours and previous researches 

on waves in such medium (e.g., Balsara 1996) is that we treat not only the propagation of waves, but also their 

generation from a source. Our formula could be extended to deal with problems concerning the global evolution 

of molecular clouds. For the first step, we present in this paper an exact solution of our formula in the case that 

the source is modeled by an oscillating source. We then apply it to suggest the role of outflows as a heating source 

and an energy source of turbulence in molecular clouds. In the next section we derive the basic equations and 

inhomogeneous wave equation. In section 3 we estimate the energy of the waves from a magnetized turbulent source. 

Applications to molecular clouds are presented in section 4, and finally we summarize the results in section 5. 

2. Basic Equations 

We consider the following two-fluid equations for a partially ionized, magnetized, and compressible medium: the 
continuity equation, the momentum equation, and the equation of state for the neutrals, 

|! + v-0H = o, (i) 

dv 

p—+p(v-V)v = -Vp + F Mc , (2) 
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and 

P = c 2 sP ; (3) 
the momentum equation for the ions, 

+ Pi(vi ■ V)«i = -Vpi + — (V x B) x B - F fric , (4) 
the induction equation, 

— = V x ( Vi x B), (5) 

and the divergence free condition of the magnetic field, 

V • B = 0. (6) 

In the above equations p, v, p, B, and c s denote the density, velocity, pressure, magnetic field, and sound speed, 
respectively. Subscript i is used to represent ions. The force .Ff r i c is the frictional force given by 

F flic = 7PPi(vi-v), (7) 

which arises from collisions between the ions and the neutrals. Here, we approximate the drag coefficient by 
7 w 3.5 x 10 13 cm 3 g" 1 s^ 1 (see e.g., Shu 1983). 

We now consider cool molecular clouds where the ionization fraction is very small (the ratio of the number density 
of ions to that of neutrals, m/n ~ 1CT 8 - 1CT 6 ) and the collision time is short compared with other concerned 
timescales. In this condition, we consider only the last two terms in equation (4), which are the magnetic force and 
frictional force, since the inertia of the ions to their motions is negligible. We, therefore, use the following equation 
instead of the equation (4): 

-^-{VxB)xB^ 1 p Pi (v i -v). (8) 

We assume that the medium is initially uniform and stationary (p = p a = const., v = vo = 0, and B = B = 
(Bo, 0,0), where B n = const.), and then waves (p = pi, v = v\, and B = B\) are generated from a compact 
turbulent source. Thus, we write the density, for example, as a superimposion of the zeroth-order undisturbed, 
uniform part, p , and the first-order perturbed part, p\\ 

P = Pa + Pi- (9) 
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In addition, we assume the isothermal condition. 

Making use of the x-component of the perturbed vorticity £ and the divergence of the perturbed velocity rj given 

by 

< 10 > 

and 

= dv lx dvi y dvu , n . 
dx dy dz 

we can eliminate v\_, pi, and B\ from equations (1) to (9). We thus obtain the following linearized inhomogeneous 
wave equations with a radiation source Q^(x,t) and Q v (r,t) (see the Appendix): 



2 1 



d 2 c\ d 3 ^ & 

dt 2 ipidtdx 2 ° A dx 2 



£{x,t) = Qfat) (12) 



and 



r) 4 r 2 f) 3 r) 2 r 2 r 2 r) r) 2 1 

° C A vy2 ° i 2 , 2\v72 ° . C A C s vji ° , 2 2^2 ° 

^ ~ ^T V A73 ~ ( C A + c s)V ^2 + — — V j- + C A C S V — T 



v(r,t) = Q v (r,t), (13) 



where ca = B /(4npo) 1 ^ 2 is the Alfven speed. We note that the wave equation for £ = dv\ x /dx is written in the 
same form as equation (13) (see equation (A18) in the Appendix). Equation (12) represents the wave equation for 
the Alfven waves and both equations (13) and (A18) are for fast and slow waves. From these three equations, we 
can derive V\ x ,vi y , and v\ z (Lighthill 1960). The source terms of equations (12) and (13) are derived as 

Q e (x,t) = [V x q] x (14) 

and 

Q,(r,i) = |kv-<Z)-C A V 2 |^ (15) 
where 

d 

q(r, t) = —g t [{vt ■ V)t> t - (V x h t ) x h t ] - c A {V x [V x (v t x h t )}} x e x . (16) 

Here, v t and h t = B t j '(47rpo) 1 ^ 2 denote the velocity and the Alfven velocity in a turbulent source, respectively, and 
e x = (1,0, 0) represents the direction of the uniform magnetic fields. We assume that the gas is fully ionized in the 
turbulent region in deriving the above source terms (e.g., in outflow lobes the temperature of the gas is high enough 
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to be fully ionized). We note that q, which comes from the non-linear terms in the basic equations, represents the 

fluctuation of the momentum flow in the turbulent source (see Kulsrud 1955; Kato 1968). In the case of a fully 

ionized plasma (where Ff r [ c = 0), the wave equations become 



dt 2 ° A dx 



£(x,t) = Qfat) (17) 
and 

- r^4 pp. pp. ' 

v(r,t) = Qr,(r,t) (18) 



d i d 2 d 2 ' 

--(d + cD^^+clclV 2 -^ 



dt 1 s/ dt 2 A 3 dx 2 _ 

(see Kulsrud 1955; Lighthill 1960: Kato 1968). 

We now assume that the oscillating source has a fixed frequency ojo; that is, the source terms are expressed as 

Q^(x,t) — Qz(x) exp(— iivot) and Q v (r,t) = Q n (r) exp(— iw t). In this case, from equations (12) and (13), £ and r\ 
are written as follows by means of Fourier transformation: 

foe-*** r Q^)eMik^) dkx (1Q) 

and 

V(r,t) = v0 e-^ f te^) d 3 fc| (20) 
Joo D v (k,ui ) 

where 

D ( (k x ,u>) = uj 2 + i-^-ujkl -k 2 x (21) 

and 

D v (k, u;)=cv 4 + *^ 3 k 2 - (c| + c 2 )uj 2 k 2 - i^ook 2 + c 2 A c 2 s k 2 k 2 x , (22) 

and Q^(k x ) and Q. q {k) are the Fourier transformations of Q^(x) and Q v (r), respectively. Here, D^(k X7 w) = 
represents the dispersion relation for the Alfven waves and D rj (k, w) = is for the fast and slow waves (cf. Balsara 
1996). Integrating equation (19) by the residual method, we can obtain £ as 

«*,*) = £oe— • 2™ V Res giMg^M , ( 23 ) 

where k xp represent poles of the function in the integral (19), which satisfy D^(k xp ,Wo) = 0. Now, only the 
wavenumber which satisfies $t(k xp ) > is used, since we consider that the waves are going forward. We can obtain 
r\ in the same manner (see Morse, Feshbach 1953; Lighthill 1978). 
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In figure 1, we show the real part of (a) £(x) = e tuJot £,(x, t) and (b) rj(r) = e lulat r](r, t) in the case that Q^(k x ) and 

Q n (k) are constant [which means that Q^(x) and Q v (x) are Dirac delta functions, representing the limiting case 

that the size of the source I — > 0]. We calculated rj{r) in the plane x = 0. In the calculations we put c s = 0.2 km 

s _1 and ca = 0.45 km s _1 (which corresponds to B = 10 pG and n = 10 3 cm -3 ). £(x) and rj(r) were normalized 

by £(0) and 7/(0), respectively, and x and r were in the units of pc. The dashed lines, solid lines, and dotted lines 

represent the waves from the source with frequencies of ui /lPi — 0.5, 1.0, and 5.0, respectively. In order to obtain 

k p satisfying D n (k p ,uJo) — 0, we used Mathematica. We note that the spatial distribution of the amplitudes of the 

waves do not vary with time, since we now consider that the source oscillates permanently with a certain frequency. 

Then, figure 1 indicates that when the frequency of the source is low, the waves dissipate less. In contrast, the 

waves dissipate soon if the frequency is higher than the collision frequency, 7p; (see also the discussions for figure 3 

in section 3). 

3. Energy of Generated Hydromagnetic Waves 

In this section we estimate the energy of hydromagnetic waves generated from a compact turbulent source. The 
energy of the waves is estimated as 

>2 / „ \ 2l 



^wave ^ 



V 



47T/90 \P0 



dV, (24) 



where vf = v 2 x +v 2 y +vf z and B\ — B 2 x +Bf y +B 2 z . Because the dependence of the variables on the distance from the 
source is the same as 77, we approximately calculated the energy as E oc J rfdV. We put dV w A-Kr 2 dr ~ Airc\t 2 dt, 
since the head of the wave propagates with about Alfven speed (which means r ~ cpj:)- Figure 2 shows the time 
evolution of the energy of waves generated from an oscillating source with a frequency uj /lPi = 10~ 8 (bold solid 
line), 0.5 (dashed line), 1.0 (solid line), 2.0 (dot-dashed line), and 5.0(dotted line). The normalization of the energy, 
Eq, is the energy provided to the waves by the source during 10 5 years. Figure 2 suggests that the energy of the 
waves is proportional to the time if the frequency of the source is very small. On the other hand, when the frequency 
is higher than the collision frequency of an ion with neutrals, the time evolution of the energy has an upper limit. 

These results can be understood from the dispersion relations of the hydromagnetic waves. In figure 3, we show 
the imaginary part of the wavenumber, hi, as a function of the frequency of the source, wq. The solid line, dashed 
line and dotted line denote the fcis of the Alfven waves, the fast waves, and the slow waves, respectively. When the 
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frequency approaches zero, k\ also becomes zero for all three waves. Since the waves are proportional to exp(ik x x) 
or exp(ifcr), as can be seen in equation (23), k\ can represent the dissipation rate of the waves. Thus, in the limit 
of vanishing frequency, there is no dissipation, as shown in figure 3. This leads to the strong conclusion that if the 
mechanical luminosity of a source has a constant value during its duration, the generated energy of the waves is 
nearly conserved; then, the generated total energy is proportional to that duration. On the other hand, k\ becomes 
larger as the frequency increases, which means that the waves soon dissipate. Then, the energy scarcely increases 
over the time when the head of the waves reach a certain radius from the source. Thus, the upper limit appears as 
seen in figure 2. 

In summary, the ratio of the dissipated energy to the energy of the propagating waves varys depending on the 
frequency of the source, and evolve with the time. The ratio (dissipated/propagating) becomes larger as the frequency 
becomes smaller or as time passes. 

4. Discussions 

In the previous section, we estimate the energy of the waves generated from a compact turbulent source. Here, 
we apply the results to the case in molecular clouds, where the outflow lobe associated with a young stellar object 
is a good candidate for a compact turbulent source (e.g., Nomura et al. 1998). 

4-1. Energy of Waves Generated from Outflows 

We first calculate the characteristic frequency and the energy of the outflow lobes in a molecular cloud, u>fl ow and 
-Eflow as 



where itfl ow , la ow , and ifl ow are the typical velocity, size, and lifetime of the outflow lobe, respectively, and Lq ow 
represents the typical mechanical luminosity of the outflow lobes in a molecular cloud. When we take the collision 
frequency, 




(25) 



and 




(26) 



No. 3] Evolution of Disturbances in Plasma 9 

as the normalization of the frequency, the frequency of the outflow lobes can written as Wfl ow ~ 2 7/^. 

Next, we estimate the energy of waves generated from the outflows. Kato (1968) studied the generation of Alfven 

waves from a turbulent source in a fully ionized, magnetized plasma. He showed that the efficiency of the energy 

inputs from the source to the waves is almost unity if the Alfven Mach number of the source is larger than unity. 

Since in the outflow lobes the velocity (~ 10 km s _1 ) seems to be super Alfvenic, its entire energy will turn into the 

waves. Thus, E used in section 3 corresponds to Efi OV! in this case. The ratio E wave /-Eflow at t = 10 5 yr with the 

source of the frequency of 2 -ypi is calculated as ~ 0.73 (see also figure 2). Therefore, the energy, 

£wavc ~ 0.73 Efi ow ~ 3.4 x 10 47 erg, (28) 

propagates as waves, and the remaining fraction, 

Efl„w - £wavc ~ 1.2 x 10 47 erg, (29) 

dissipates and heats the neutrals due to ion-neutral damping. 

At last we compare the timescale of the energy injection from the outflows, t; n j, with the timescale of the wave 
dissipation via the ion-neutral damping, ij-jj;™. Because the injection timescale is comparable to the lifetime of an 
outflow, fflowi we obtain 

i inj ~ 1 x 10 5 yr. (30) 
Because the dissipation timescale corresponds to the collision time, (7Pi) _1 , 

t^~2xlO«yrf^yY- J A =ff y 1 . (31) 



W~' J V 10 cm- 
Prom equations (30) and (31), we find that i; n j is slightly shorter than tj^ e ; that is, ii„j ^,t™ s c . This means just 
the same as the above result; i.e., more than a half of the source energy propagates as waves. 

In the following, we compare the above-mentioned energy and timescales with others in molecular clouds to 
consider a significance of outflows as an energetic source in molecular clouds. 

4-2. Energetic Source of Interstellar Turbulence 

Here, we consider the energy inputs from outflows to turbulence in a molecular cloud. The absorption or scat- 
tering of waves propagating in a non-magnetized turbulent medium have been investigated both theoretically and 
experimentally (Hunter, Lowson 1974; Noir, George 1978; Salikuddin et al. 1988). In a partially ionized, magnetized 
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medium, hydromagnetic waves may be absorbed in the same manner. We thus compare the energy of the outflows 

with that of the turbulence in a molecular cloud to examine the possibility for the activity of the outflows supporting 

the turbulence. The energy of turbulence in a molecular cloud, £Wb is calculated as 

where M c i ou d and cr tU rb are the typical mass of a molecular cloud and the turbulent velocity at the typical size of 
the molecular cloud, respectively. From this and equation (28), we find that the energy of interstellar turbulence is 
comparable to the energy of hydromagnetic waves in a molecular cloud, namely, -EWrb ~ E wave . 

By recent numerical simulations, the dissipation timescale of interstellar turbulence were calculated as 

i^ r s b ~ (0.25-2) x 10 6 yr (33) 

(Stone et al. 1998). Equations (30) and (33) suggest ti n j<,t^l^, which means that the energy inputs from the 
outflows to the interstellar turbulence can balance with the dissipation of the turbulence, which occurs mainly due 
to the shocks. Thus, the turbulence may be in equilibrium in molecular clouds. In any case, the outflows play an 
important role for energy inputs to the turbulence in molecular clouds. 

4-3. Heating Source in Molecular Clouds 

Here, we confirm that the dissipation of waves generated from outflows plays an important role in molecular 
clouds. A variety of heating sources have been considered in molecular clouds, such as cosmic rays, H 2 formation on 
grains, gas-grain collisions, the compression or collapse of dense cores, and dissipation of turbulence, etc. (Spitzer 
1978; Hollcnbach 1988). The heating rate due to the cosmic rays, which is common as a heating source in molecular 
clouds, is 

r CR ~ 6.4 x 10- 25 crg cm-V 1 (^-^—^J . (34) 

Stone, Ostriker, and Gammie (1998) have estimated the heating rate based on the dissipation of interstellar turbu- 
lence, making use of numerical simulations, as 

^turb \ / ^cloud 
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where l c ioud is the typical size of a molecular cloud. We now calculate the heating rate via the dissipation of the 

waves generated by the outflows as 



-i -, , , : / ^flow -E'wavo \ / ^flow \ ( ^cloud \ ' lna\ 

1.5 x 10 erg cm °s j= — = — — . (36) 

M.2 x 10 47 erg/ V10 5 yrj V 10 P C 



This is slightly larger than the heating rate of the other heating sources. We note that this heating rate is comparable 
to the cooling rate, due to optically thin CO line emission, 

A cool ~ 4.1 x l(T 24 crg cm-V 1 (y^) ' ( 37 ) 

where T is the temperature in the molecular clouds (Goldsmith, Langcr 1978). We thus conclude that dissipated 
waves originating from outflows is one of the heating sources in molecular clouds. 

5. Summary 

We investigated the evolution of hydromagnetic waves generated by a compact turbulent source in partially ionized, 
magnetized plasmas, by applying the Lighthill theory. We calculated how much of the energy of the generated waves 
is dissipated to the ambient molecular gas via ion-neutral damping, depending on the frequency of the oscillating 
source. As a result, we found that the ratio of the dissipated energy to the energy of the propagating waves becomes 
larger as the frequency becomes higher, and also increases with passing time. The new point of our paper that is 
superior to Balsara (1996) is that we treat the wave equations with source terms as Lighthill did; namely, we consider 
the generation of waves from a source as well as their propagation. 

The implications to the phenomena in molecular clouds are also presented. We discuss the effects of outflows as 
an energy source in molecular clouds when they act as compact turbulent sources. Our conclusions are: 

1. Since the timescale of the ion-neutral damping is slightly longer than the typical lifetime of the outflows, more 
than half of the energy generated by the outflows propagates as waves, and the remainder dissipates in the molecular 
gas via ion-neutral damping. 

2. The energy of the propagating waves is comparable to that of the turbulence in molecular clouds. Interestingly, 
the dissipation timescale of the turbulence is comparable to the lifetime of the outflows, suggesting that it is possible 
for the hydromagnetic waves generated by the outflows to support the turbulence in molecular clouds. Although 
our discussion partly confirms the suggestions by Myers and Lazarian (1998), the essential difference between ours 
and the previous studies is that we consider the evolution of waves originating from the outflows. 
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3. The heating rate due to outflows through the ion-neutral damping of the generated waves is slightly larger 

than those of the other well-known heating mechanisms in molecular clouds. That is to say outflow is one of the 

dominant heating sources which sustain the temperature of the molecular clouds. 

To conclude, we confirm that outflows generally play a substantial role as an energy source in molecular clouds, 

in the context of the Lighthill theory Although we have applied our formula to only simple cases for the first step 

in this paper, it could be developed to contribute to studies on the global evolution of molecular clouds. 

We would like to thank the referee, Prof. S.W. Stahler, for his careful reading of this manuscript and valuable 
comments, which improved the clarity of our discussions. 



Appendix. Derivation of Equations (12) and (13) 

Making use of equation (9), the basic equations are written as follows: the continuity equation, 
dpi ( dv lx dviy dvi z \ 

ln +p \-e^ + l£ + li7)=^ (A1) 

the momentum equation with equations (3), (7), and (8), 



+ [(vi ■ V)vi] x = -f^ + [(V x h) x h] x , 
dviy r , . . c 2 s dpi B Q f dB X y dB lx 



(A2) 



dt + ^ ' V ^ ^-T^l^f if) + r * ^ x H, 

and 

dvi z ,, „. . ct dpi B ( dBi v dBi x , . ._. 

+ [(vi ■ V) Vl ] z = ~ — + -j— ^- - )+[(Vxh)x h] z , (A4) 



dt po dz 47TP0 V ® x 

where h = Bi/^irpo) 1 / 2 , and the induction equation with equations (6) and (8), 

dBi x _ Bo (^y + ^A +[Vx{BiX Vl)]x = A^B 1X , (A5) 



dt \ dy dz J -fp 

dBi y 



dt 

and 



„ dviy r „ . c A d ( dBi v dBi x \ , , . 

- B ^ + [V x (Bi x Vl)]y = + m - j , (A6) 



dB lz dv lz c\ d (dB Xz dB lx , 

~dT - + [V x (Bi x vi)) z = —-{—- — ). (A7) 
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At first, we derive equation (12). From equations (A3) and (A4), and the variable, 

_ dv lz 8v ly 



dy 8z 

we obtain 



d 2 Z 3 fr7 ., ^ N n B Q 8 ( d 8B ly d dB lz \ 8 r „ u ^ ,_, Lll 



(A8) 



(A9) 



By substituting equations (A6) and (A7), and using the assumption that the gas is fully ionized in the turbulent 
source region, the first term on the right-hand side of equation (A9) becomes 



Thus, we can derive equation (12), 



8t 2 iPidtdx 2 ° A dx 



£(x,t) = [Vxq} x , (All) 



where 



8 

q{r, t) = -Q- t i(vi ■ V)«i -(Vxh)xkj- c A {V x [V x ( Vl x hi)]} x e x . (A12) 

Here, e x = (1, 0, 0) is the direction of the uniform magnetic fields. 

Next, we derive equation (13). From equations (A1)-(A4), and the variable 

v= d^ + d^ + ^ 
ox ay 8z 

we obtain with equation (6), 

g + |{V ■ • V). l]} = c 2 ^ + ^-^B lx + -§{V • [(V x h) x h]} . (A14) 
By substituting equation (A5), the second term in the right hand side of equation (A14) becomes 

4V 2 (?7 - C) + — t^V 2 ^ - c A V 2 [V x(hx vi)] x , (A15) 

where £ = dv\ x /dx. Now, from equations (Al) and (A2) we obtain 
rP-C 8 8 f) 2 n 8 8 

W + SS><- - VKI " " c - a? - 5 &1< V * »> * *!- < A16 » 

Thus, making use of equations (A14)-(A16) and the fully ionized condition in the source region, we can derive 
equation (13) and wave equation for (, 



14 
and 
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«4 2 Q3 x2 2 2 o a2 1 

dt± ' 7 Pi dt* (CA + Cs)V 9i 2 + 7 Pi ^ + CaCs V fe 2 



((r,t) = Q c (r,t), 



(A18) 



where 



<9 



Qr,(^*) = ^2(V-q)-4V 



2 rj2 



3a; 



(A19) 



and 



r «2 
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Figure Captions 

Fig. 1. Real parts of (a) the x-component of perturbed vorticity, £(x), and (b) the divergence of perturbed velocity, 
rj(r), generated by a compact source with frequency of ujo/jpi — 0.5 (dashed lines), 1.0 (solid lines), and 5.0 
(dotted lines). 7/5; represents the collision frequency, and the compact source is located at x = 0. 

Fig. 2. Time evolution of the wave energy generated by a source with a frequency of wo/7Pi = 10~ 8 (bold solid 
line), 0.5 (dashed line), 1.0 (solid line), 2.0 (dot-dashed line), and 5.0 (dotted lines). Clearly, when the oscillation 
frequency is larger than the friction frequency, 7/5;, the wave energy dissipates soon in the molecular gas via 
ion-neutral damping. 

Fig. 3. Imaginary parts of the wavenumbcr, fci, as a function of the source frequency, uiq- The units of the frequency 
is friction frequency, 7pj. fci for three types of the hydromagnetic waves, that is, the Alfven mode (solid lines), 
the fast mode (dashed lines) , and the slow mode (dotted lines) are represented. When the frequency of the source 
is large, the generated waves easily damp, because their amplitudes are proportional to exp(— |fci|x). 
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